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Abstract 

In this note we show that the random homological product code con¬ 
struction of Bravyi and Hastings can be extended to qudits of dimension 
D with D an odd prime. While the result is not surprising, the proof does 
require new ideas. 


1 Introduction 

Recently it has been shown that the homological product of two random chain 
complexes gives rise to good CSS codes with stabilizer weight 0{^Jn), where 
n is the code size. mm The proof makes use of a simplified version of the 
homological product, called the “single sector theory”. In this construction 
a chain complex is a pair (C, d) with C a free module over a ring R with a 
distinguished basis and the linear map d : C —>■ G satisfies = 0 . When 

2 = 0 in i?, it is possible to define tensor products of such chain complexes as 
{Cl ® C2,di® I2-\- Ii® d2)- It is then shown that if i? = Z2 and di and 82 are 
random conjugates of some fixed boundary operator, then the product complex 
leads to good codes with high probability. 

Unfortunately, when 2 ^ 0 , di ® I2 + Ii ® 82 is not a boundary operator 
in general, since its square is 2 di ® 82- In typical applications of homological 
algebra, the modules Ci and C2 are graded and boundary maps are homoge¬ 
neous of degree —I (or -|- 1 ), and so it is possible to remedy the situation by 
introducing signs depending on the grading. In ref. [ 2 ] this variant is called the 
“multiple sector theory”. However, in this case it seems to be difficult to find 
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good lower bounds on the distance of the homological product code. Our main 
result is the analysis of an intermediate structure, one which is rich enough to 
make it possible to form tensor products even in the qudit case, and at the same 
time it is sufficiently close to the single sector variant so that the proofs can be 
modified to reach essentially the same conclusion. 

The structure of this paper is as follows. In section [ 5 ] we describe the above 
mentioned intermediate structure. In section | 3 ] we extend the proof of Bravyi 
and Hastings to qudit codes. Here we mainly focus on proofs that are signifi¬ 
cantly different from those of the single sector theory. While the main line of 
argument is presented here, the reader may refer to [ 5 ] for further details. 


2 Double sector theory 

To define the “single sector” theory, one starts with a single Z^i-module C with 
a basis and a boundary operator S : C ^ C, which satisfies 6“^ = 0 . In this 
case the columns and rows of (the matrix of) S give Z and X-type stabilizer 
generators. 

When D = 2 , the tensor product of two chain complexes {Ci,Si) can be 
defined in the single sector theory simply as (C, d) with C = Ci ^ C2 and 
d = Si (Si I2 + h (S' 62- This actually works over any ring as long as 2 = 0 . Ref. 
[1] uses the single sector theory to form the product of two independent random 
chain complexes, because in this framework the distance of the resulting CSS 
code is easier to analyze. 

When U > 2 , it is not possible to form a tensor product of chain complexes 
in the single sector variant. However, it is possible to introduce some extra 
structure which lies between the single and multiple sector theories in the sense 
that it allows us to form tensor products, but it is not much more difficult to 
find lower bounds on the distance of the product of suitably defined random 
codes. One possibility is given by the following definition (see e.g. 0 ): 

Definition. A differential Z^i-module with involution (or chain complex with 
involution over is a triple {C, 9 , P) where C is a module over d : C ^ C 
and P : C ^ C are Z^j-linear and satisfy 

9 ^ =0 

P^=I (I) 

dP + Pd = 0 . 

The tensor product of two differential Z^i-modules with involution (Ci, 61, Pi) 
and {C2, S2, P2) is defined to be (C, 9 , P) where C = C'i(8)C'2, 9 = 5i®l2-\-Pi®S2 
and P = Pi® P2. 

With this definition, the tensor product also satisfies eqn. ([T|). When D = 2, 
the single sector theory suffices, whereas when D is not a prime, additional 
difficulties arise, and we do not know how these can be overcome. For this 
reason, from now on we will assume that D is an odd prime. In this case Z^i is 
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a field, therefore the appearing modules are in fact vector spaces. In addition, 
2 is invertible, so we can form the linear combinations 


P+ = 


I + P 
2 


and P- = 


I - P 
2 
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where / is the identity operator. These satisfy P^ = P+, P^ = P- and P+P- = 
P-P+ = 0. Moreover, P+ + P_ = / and P+ — P- = P. This implies that C can 
be written as a direct sum C = C+ © C- with C± = P±C, these are invariant 
subspaces of P and P\c± = ±idc±. It also follows that dP± = P+d, therefore 
d takes C± to C+. Equivalently, we could have started with a pair of spaces 

C+ and C- and a pair of linear maps _: C_ —>■ C+ and d _: C+ —>■ C_ 

satisfying _ d _= d _|-9-|_= 0. In line with the naming convention of ref. 

[5] , this may be called the “double sector theory”. 

The key feature of the homological product is that it does not increase the 
stabilizer weights too much. Let w^A) denote the maximal weight over the 
rows and columns of a matrix A (when A is a linear map, this depends on the 
basis, but we will always use a fixed basis). Clearly, w{A + B) < w{A) + w{B) 
and w{A ® B) = w{A)w{B) holds with respect to the product basis. In the 
single sector theory these imply that w{d) < w{Si) + w{S 2 ), but in our case this 
relation becomes 

w{d) < w{5i) + w{Pi)w{S 2 )- (3) 


In order to get the best possible bound, we will choose Pi to be (in block matrix 
form) 


P^ = 


I 

0 


0 

-I 


(4) 


Then w{Pi) = I, therefore the same inequality holds as in the single sector 
theory. An additional advantage of this choice is that C+ and C- are now 
spanned by subsets of the basis vectors. In the next section we would like to 
consider multiple boundary operators on a single space. One possibility is to 
choose a hxed boundary 5 and conjugate it with arbitrary invertible matrices. 
In order to retain compatibility with the involution, we need to conjugate that 
with the same matrix as well, but this would potentially increase its row and 
column weights. This can be prevented by requiring that the invertible matrix 
commutes with P. In a block matrix form, the commutator reads 


I 0 
0 -/ 

Since 2 is invertible, the commutator is 0 iff P = 0 and (7 = 0. The conjugating 
matrix is therefore block diagonal. Let us call its diagonal blocks U+ and C/_. 
Then the boundary operator transforms as 


'U+ 

0 ■ 


■ 0 

<5+-' 


'U-^ 0 


0 U+5+-UZ^' 

0 

U-_ 


+ 

1 

0 


0 c/r\ 


U-5-+U+^ 0 



A 

B 


A 

B 


I 0 


0 

2B 


c 

D 


C 

D 


0 -I 


-2(7 

0 


It is possible to show that both d-i_and 5 _^ can be brought to a standard 

form simultaneously, which only depend on the two homological dimensions. 
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We will choose random boundary operators by fixing a standard do and letting 

d-i_= U+SoUZ^ and S _|- = U-SqU^^ with {7+ and U- drawn uniformly from 

the set of invertible matrices. We choose do to be 


do = 


0 

0 

0 


0 

0 

0 


0 

I 

0 


(7) 


in block matrix form with row/column sizes H, L, L in this order, where H is 
the homological dimension and H + 2L = dimC_|_ = dimC_. In particular, we 
choose these parameters to be the same for the + and — blocks, but it would 
be possible to allow different values. 

Once we have two chain complexes with involution and form the tensor 
product, it would be possible to forget the involution and construct a CSS code 
as in the single sector theory. Unfortunately, the code distance of the product 
of two random complexes seems to be difficult to bound in this case. Instead of 
this, we adopt a definition which resembles the multiple sector theory. Phyisical 

qudits will correspond to the basis elements of C+, the columns of _will be 

Z-type generators and the rows of d _^ will be the X-type generators (for this 

we need the special form of P given above). The code obtained this way will be 
denoted by CSS{C,d,P). 


3 Homological product of two random chain com¬ 
plexes with involutions 


In this section we let n = + 2L with P[ = \_pn\ (where p is related to the 

coding rate), and So is the matrix from eq. ([7]). By a random 2n x 2n boundary 
operator on C = C+ 0 C- with C± = we mean 


'U+ 

0 ' 

'0 

<5o' 

0 

0 

U-_ 

.<5o 

0 

0 uz\ 


in block form. The boundary operator in the middle will be called the standard 
one. 


Lemma 1. For any e > 0 there exist c, p > 0 such that the prohahility that the 
kernel of a random 2n x 2n boundary operator with H < pn contains a nonzero 
vector with weight less than cn is 0 (ZI*'“^ 72 +e)n^^ 

Proof. Let S be as in eq. (|8]). Since 


U = 


U+ 

0 


0 

U- 


(9) 


is invertible, kerd = Ukei{U ^SU). P anticommutes with 6 , therefore if u G 
kerd then v± = ^{I±P)v G kerJ, which in turn holds iff U~^v± G ker(P~^SP). 
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For any nonzero vector w € ker([/“^i5?7) there are three possibilities accord¬ 
ing to whether the vectors w± = ^(J ± P)w vanish or not. If e.g. ^ 0 and 
W- = 0, then Uw is distributed uniformly on the set of vectors v with 7 ^ 0 
and = 0. When c < 1 — the probability that such a vector has weight 
less than cn is at upper bounded by 

^ 23(-(l-‘=)'°gi3(l-c)-'=logo c-|-clogD(D-l)-l)ra^ 

( 10 ) 

We get the same probability when w+ = 0 and W- 7 ^ 0. When w+ 7 ^ 0 7 ^ W-, 
the vector Uw is uniform on vectors v with v± ^ 0. This time the probability 
of having weight less than cn is upper bounded by 

0(1) • logi3(l-c)-clogo c+cloBo{D-l)-l) 2 n ^ 



Finally, we apply the union bound. The number of vectors of the first two 
types is 0 ( 1 ) • 02"+2 Lp "1 and the number of vectors of the third type is 0 ( 1 ) • 
jjn+ lpn\ ^ ^];j 02 . 0 fore the probability of having a nonzero vector with weight less 
than cn is at most 

0 ( 1 ) • c-|-clogo(D-l)-i-|-f )ra_ 

For any O > 0,e > 0 and small enough c,p the exponent is less than (—5 + 
e)n. □ 

Lemma 2. Let {C,d^P) be a chain complex with involution, V <C a subspace 
and a direct complement. Let W, : C ^ C be the projections corre¬ 
sponding to this direct sum decomposition and suppose that PW = WP. Let 
= W5(y^) and V' = V /. Let ip : C ^ V' be defined as h ^ Wh-\- . 

Then the maps 6 ', P' : V ^ V' given by 

5'{x -\- S^) = (p{d{x)) and P'{x -\- S^) = (p{P{x)) (13) 

are well defined, (V^ 6 ', P') is a chain complex with involution, and ip : C ^ V' 
is a chain map satisfying ipP = P'ip. Moreover, the equalities 

kerJ' = i^(J“^(V^)) and imd' = (/?(im^) (14) 


hold. 

Proof. Suppose that x G , i.e. there is some y G such that x = WSy. 
Then WSx = WSWSy = W5{W — I)5y = —WSW^Sy S 5'^, therefore ip{5x) = 
W6x . Similarly, ip{Px) = tp{PW5y) = tp{—W5Py) G , because 

Py = PW>y = W>Py G V>. 

By definition, S'(p = (p6 and P'ip = ipP. This implies that 5'5'ip = ip55 = 0, 
P'P'if = ipPP = ip and {P'6' -G 5'P')ip = ip{P5 -f 5P) = 0. By surjectivitiy of 
ip, the equalities 5'^ = 0, P'^ = I and P'5' -G 6'P' — 0 follow. 

The proof of eq. m is the same as in Lemma 8] (Lemma 10 in [T]). □ 
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Definition. A boundary operator S : C+ 0 C- —>■ C+ 0 C- is called good if 
neither ker 6 fl (7+ nor ker 6 fl C- contains a nonzero vector supported on the 
last n — n' coordinates. 

In the following V will be the subspace of C spanned by the first n' basis 
vectors in C± and the direct complement spanned by the remainig ones. 
Then the conditions of lemma [5] are satisfied. 

Lemma 3. Let 6 be a good boundary operator. Then dim = 2n' — n and 

dim(ker(5^_i_) = dim(keri5 |_) — (n — n') (15) 

dim(im(5(__|_) = dim(im(5 |_) — (n — n') (16) 

hold, and similarly for (5(_. 

Proof. The proof is very similar to that of [U Lemma 9] (Lemma 11 in [T]). □ 

Let Cl = C 2 = equipped with Pi = P 2 which is diagonal in the standard 
basis, and acts as 0l (—1) on the first (last) n coordinates, and let 5 i ,52 be 
compatible boundary operators. Suppose that if G (7+ fl ker 9 where C = 
Cl 0 C 2 , C+ = Ci+ 0 C 2 + 0 Cl- 0 C 2 - < C and 9 = di 0 /i 0 Pi 0 <52- "0 can 
be thought of as a block-diagonal matrix with two n x n blocks V'-i- and if-. 

If the weight of if is less than cn^ then the same is true for both blocks. 
Choosing some r such that c < r < 1, we can find at least n' = (1 — r)n rows 
and columns in if± having weight at most cnr~^. The two n' x n' matrices 
obtained this way will be called the reduced matrix of if. An n' x n' matrix 
having row and column weights at most c'n' (with c' = cr“^/(l — r)) is said to 
satisfy the uniform low weight condition. 

In the following lemma, Vi and Vf" are the subspaces of Ci defined similarly 
as before. 

Lemma 4. Suppose that the distance o/CSS((7i, Si, ±Pi) is at least 2{n — n') +1 
(i = \,2). If h G (7+ nker9 has vanishing reduced matrix then /i G im9. A 
similar statement holds for cocycles. 

Proof. The proof goes along the lines of [2l Lemma 5] (Lemma 5 of [I]), but 
there are differences due to the fact that we restrict to the subspace (7+. Let 

ha+ G ker(5j__|_ \im5j^_and ha- G ker5j_^ \im5j^_be nontrivial cocycles 

(a = 1, 2) and let S = {n' + l,n' + 2,... ,n} U {n + n' + l,n + n' + 2,..., 2n}. 
Since [S’] = 2(n — n') < 2{n — n') + 1, the cleaning lemma ([H Lemma 1], 

the same proof works for qudits) implies that there exist uJa+ G imJj,^_and 

uia- G im(5j__|_ such that the support of ha± + LOa± is disjoint from S. 

Choose nontrivial cocycles supported outside S such that their cosets span 
the cohomology groups: 

ker5j_+ =span(hi+,/i2^,...,/if+) 0im(5j+_ 
ker = span(/ii _ ,hl_,...,h^_)+ im . 
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Take the dual basis of nontrivial cycles: 

ker(5a-+ = span(/i^+,...,+im(5a+- 
ker 5a+- = span(/i^_, ,..., /if_) + im 5^-+, 


(18) 


such that (/ia±, h^a±) — ^ij- By the Kiinneth formula we have 

ker9 n C+ =span{/i^_|_ 0 <i,j< H} 

+ span{/i 2 - ® < b j ^ H} + imi9 n C+ 

keid^ n C+ =span{/ij_|_ 0 < i,j < H} 

+ span{h^_ 0 h 2 _|l < *, j < H} + imi9^ fl C+. 

Let h G ker d fl C+ be a cycle with vanishing reduced matrix. Then 

H H 

h = ^ 0 ^ x~jh\_ 0 h{_ + w, (20) 

*j=i *j=i 

where w S im9 fl C+ and = (h\± 0 (i 2 ±,/i) G Z^. Since h\^ 0 is 


supported on the reduced matrix where h vanishes, xf^ = 0, therefore h = uj G 
imd. 

The proof for cocycles is the same with the roles of S and 6 ’^ reversed. □ 

Definition. We let denote the number of rank R matrices of size Ax B 

which extend an (arbitrary) rank r matrix of size axh over Zd. 

Lemma 5. 

^a,b,r = 0 ( 1 ) . J^{A+B)R-R^ ^ ( 21 ) 

and 

J^A,B,R _ Q( 1 ) . £)iA+B-b)R-br-R'^ + (b-a+r+R)^/4 (22) 

Proof. The proof is very similar to that in Appendix A] (Proposition 1 in 

m- □ 


Lemma 6. Let Si, 62 be boundary operators with dimimSa = 2 L and dim ker = 
2{L + H) and let 9 = 0 /2 + Pi 0 (52- Let ZH,L{x+,r-) be the number of 

h G ker 9 ("1(7+ such that rkh± = r±. Then ZH,L{r+,r-) only depends on 
r+,r-,H and L and 

2L 

ZhA^+A-) < 0(n) ■ i:)2(ff+i)0++r-)-(r^+P ) j^-l^ + (r++r—2H)l _ (33) 

1=0 

Proof. Since di and S 2 are conjugates of the standard boundary operator 6 and 
left and right multiplication by a block diagonal matrices Ui,U 2 does not change 
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the ranks, we may assume that 5a are the standard ones. In this case (7+ Piker9 
is equal to the set of matrices of the following form: 



1_1 



B+ 

O' 


'A_ 

B_ 

O' 

h = 

where = 


D+ 

F 

,h- = 

C_ 

D_ 

G 



0 

G 

0 


0 

-F 

0 


To count the number of such matrices with given pair of ranks, let us first 
fix F and G with / = rkf and g = rkG. These ranks must satisfy 0 < 
f + g < min{r+,r_} and f,g<L, but are otherwise completely arbitrary. Let 
U, V, X, Y be invertible Lx L matrices such that UFY and VGX have nonzero 
elements only in their upper left / x / and g x g corners, respectively. Apply 
the transformations 



■/ 

0 

o' 


■/ 

0 

O' 


'/ 

0 

o' 


■/ 

0 

O' 

h-\- 1 — y 

0 

u 

0 

h+ 

0 

A 

0 

, h— 1 — y 

0 

V 

0 

h- 

0 

Y 

0 


0 

0 

V 


0 

0 

Y 


0 

0 

u 


0 

0 

A 


(25) 

This does not change the block structure or the ranks. Removing the rows and 
columns from h+ and h- corresponding to the nonempty rows and columns of 
the new F and G decreases both ranks by f+g independently of A ±, B ±, C ±, D±. 
The remaining nonzero matrices have sizes {H + L — f) x (H + L — g) and 
{H + L — g) X {FI + L — /). These can be chosen arbitrarily as long as their 
ranks are r± — f — g. Summing over the possible choices gives 


Z//,L(r+,r_) = ^ ^L,L,f^L,L,gj^2(f+g)(H+L)-2fg 


f,g=o 

/+g<min{r+,r_} 


. ^H+L-f,H+L-g,r+-f-g^H+L-g,H+L-f,r.-f-g 

= 0 ( 1 ) • Z)‘^(^+FG++r-)-rl-F_ 


(26) 


E 


D 


-U+g)^ + if+a)ir++r- -2H) 


f,g=o 

/+g<min{r+,r_} 


Introducing I = f + g, we may replace the sum over / and g with n times a sum 
over I with 0 < I < min{2L, r+, r_} to get an upper bound. □ 

Lemma 7. Let (5i , 62 be good boundary operators with Ci± = 1/j^ and homolog¬ 
ical dimension F[ -\-H. Let r(i?+,i?_) be the number of reduced cycles {h+, h-) 
with ranks rkft,-|_ = and rk/i_ = i?_. Then 

r(R+,i?_)= ^ ^ 

r_|_—0 r_—0 


where K = 2n' — n. 


(27) 



Proof. From lemma [5] we get a map (f ® Lp ■. C ^ V'l ® V '2 with C = Ci ® C 2 
and V' as in the lemma. The proof of [U Lemma 10] works in our case without 
modification, and it implies (after restricting to the subspaces C+ = 0 

C 2 + © Cl- © C 2 - and © V 2 + © V]_ © V 2 _) that 


T{R+,R-)= ^ \{g+®g-GV+\rkg± = R± and {tp^(p)g = h}\. (28) 

hGV'j^nker d' 

A similar argument as in the qubit case allows us to rewrite the above expression 


as 


r(i?+,i?_)= ^ |{/i £ n ker9'| rk/i± 

r-i-— 0 


r±}| • E 




-ppu' ,n' ,R- 
^K,K,r- 


= E E 


ZH,L-(u-n’) (?'+ , r- 


E 


n ,n' 

K,K,r- ' 


r +=0 r _=0 


(29) 

□ 


Let Zfi^ {Ui , 1 / 2 ) denote the set of reduced cycles for d = Si (Si I 2 +Pi (^62 
with block ranks i?±, where Ui are block diagonal with two n x n blocks and 


S^ = U, 


0 

<5o 



(30) 


By definition, (17i, 172)1 = r(i?+, i?_) when the Si are good. 

Lemma 8. R is possible to parameterize the sets {Ui, U 2 ) with integers 

j = 1 ,... ,r(i?+,i?_) in such a way that conditioned on < 51,^2 being good, the 
distribution of the jth reduced cycle is uniform on the set of pairs of n' x n' 
matrices of ranks (i?+, R -). 


Proof. Let Si = Vi 


0 

<5o 


5o' 

0 


^ be good boundary operators. 


subgroup GL{n',Z d)^ < GL(n,Z£))^ of matrices 


Consider the 


U = 


[/; 0 0 o' 
0/00 
0 0 u'_ 0 
000 / 


(31) 


where the block sizes are n',n — n',n',n — n'. If Ui are such matrices, let 
U[ = © U[_. Since <5^ is good and conjugation by Ui does not affect the 

last n — n' coordinates, UiSiUi is also good and Z^^ (CiVi, U 2 V 2 ) = {U[ © 
U^)ZR+,R-iVl,V2). 

Choose an arbitrary parameterization fvi,V 2 ■ [r(71+, R-)] —>■ Z}i^^}i_ (Vi, V 2 ) 
and let fjji Vi,U2V2{j) = {U'i®U!2)fvi.. y 2 {j) for all Ui,U 2 as above. Doing this for 
one representative (Vi,V 2 ) in each coset gives a parameterization for all good 
pairs which has the desired properties. □ 
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Lemma 9. For any e > 0 there is a c' > 0 such that for any 1 < R < n' the 
probability that a random n! x n' matrix Z chosen uniformly from the rank R 
matrices has row and column weights at most c' is 0(1) ■ - 2 {i-e)n r 

Proof. Similarly as [U Lemma 14] (Lemma 6 in [T]). □ 

Lemma 10. Let (5i,(52 be random boundary operators. Let denote the 

probability that there is a nonzero reduced cycle for d = 5i ® L 2 -\- Pi ® 82 
satisfying the uniform low weight condition conditioned on 5i and 62 being good. 
Then for any e > 0 and 0 < r < 1 there is a c > 0 such that 

Pred'^ < 0(n®) • £)-2(l-2»--2e+2re-p)n_ (- 32 ) 

Proof. Choose a parameterization of the reduced cycles as in lemma [S] Then the 
jth reduced cycle in (Oi, U 2 ) is distributed uniformly on the set of pairs 

of n' X n' matrices with ranks R± . If c is small enough, then lemma [HI implies 
that the probability that the jth reduced cycle satisfies the uniform low weight 
condition (with c' = cr-V(l - r)) is 0 (1) • i:)«+- 2 (i-E)"'fl+£)iJd- 2 (i- 6 )n'_R_ _ 
Summing over every reduced cycle gives 


< 0(1) ^ r(i?+,i?_)i:)^++^~-2(i-'^)"'(^++^-). (33) 

R± — l 

Next we plug in the bounds for T, Z and E from lemmas [T] |6|and[5| Then 
we extend the sums over r± to R± and exchange the order of the sums to get 


pgood - 
-^red — 


2n' — n — H 

0(n) ^ JJ-T-2HI 

1^0 


\rj^—0 R^—max{l,r^} 


jjR\-2{l-t)n' R+ jjnR+-^R\ jj{H+R+/2)r+-^r\+lr+ 


(34) 


where the square appears because the sum over (r_, i?_) gives the same factor. 

The exponent in the inner sum is quadratic, for a fixed r+ it has a minimum 
at R+ = 2n — r_|_ — 4n(r + e — re), which is larger than the upper endpoint, 
therefore the maximum is at i?+ = max{l, r+}. If r+ = 0 then we get j + n — 
2n'(l —e), while r+ > 0 gives {H + l + n — 2n'{l — e))r+. If I > —L[ — n + 2n'{l — e) 
this is increasing, therefore the maximum is at r_|_ = n', otherwise it is decreasing 
and the maximum is at r+ = 1. We get an upper bound by keeping the three 
possible maxima and counting the number of terms (0(n^)). The square of 
a sum is less than the sum of squares times the number of terms, which is 
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bounded, therefore 


pgood 

"^red 


2n' — n — H 


< O(ri^) 'y ^ ^jj^+2n-4n'(l-e) jj2{H+l+n-2n'{l-e))n' 


1=0 


_|_ £)2(H+l+n-2n'(l-e))y 


After expanding the product, each exponent is quadratic in I and they have 
maxima at 0, 0 and n' — H, respectively. Multiplying the maximum with the 
number of terms in the sum gives 


pgood ^ —2n(l — 2r — 2e+2re) _|_ p — n^{l—4r+3r^—4e+Sre—4r^e—p^) 

+ £)-2"(l-2r-2e+2r£-p) j 


where we used H = pn and n' = (1 — r)n. When n is large, the last term 
dominates. □ 


Theorem 1. For any large enough n G N there exist CSS codes with parameters 
[2n^, cn^, 2n]. 

Proof. Let (5i,(52 be random boundary operators and d = Si 'Si h + Pi 'Si S 2 as 
before. Then CSS{C 0 C, d, Pi S P 2 ) encodes 2(pn)^ qudits and its stabilizer 
weights are at most 2n. 

The probability that there is a choice of n' rows and columns in both blocks 
such that there is a nonzero reduced cycle (on these rows and columns) obeying 
the uniform low weight condition can be bounded from above as 


^bad 


+ (l_pbad) 


n' 


pgood 

-^red 


< 0 ( 1 ) + 0 {n^) 



4 

jj — 2{l—2r — 2e-\-2re—p)n 


(37) 


bv lemmas IT] and fTOl Using that 


( ) _ £)(-»'logDr-(l-r)logo(l-i-))n+o(n) /gg\ 

V(l-r)ny 

we get that for sufficiently small r and e this probability goes to 0. The same 
bound applies to cocycles since d and 9^ have the same distribution. Therefore 
when n is large enough, there is at least one i5i, ^2 such that CSS{C SC, d, Pi S 
P 2 ) has distance > cn^. □ 
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